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It is shown that quantum fluctuation theorems remain unaffected if measurements of any kind 
and number of observables are performed during the action of a force protocol. That is, although 
the backward and forward probabilities entering the fluctuation theorems are both altered by these 
measurements, their ratio remains unchanged. This observation allows to describe the measurement 
of fluxes through interfaces and, in this way, to bridge the gap between the current theory, based on 
only two measurements performed at the beginning and end of the protocol, and experiments that 
are based on continuous monitoring. 
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In the last decade fluctuation theorems have experi- 
enced a renewed and widespread interest [H-@- These 
theorems yield rigorous predictions for nonequilibrium 
processes beyond linear response theory. In particular, 
they quantify the probability of events that are forbid- 
den by the second law of thermodynamics as being ex- 
ponentially suppressed compared to the probability of a 
typical, allowed event. 

While experimental verifications of the fluctuation the- 
orems for classical systems were performed by different 
groups and for different systems, e.g., 0-0, an exper- 
iment with quantum devices was performed only very 
recently (Io| . In that experiment the flow of electrons 
through a double quantum dot placed between two leads 
with different chemical potentials was continuously mon- 
itored and the quantum fluctuation theorem was veri- 
fied for the probability p(q) that a number q of electrons 
is exchanged by the leads in a certain interval of time 
t = tj — t , see Fig. 1. With regard to theory, however, 
the existing derivations of fluctuation theorems for quan- 
tum transport do not allow for a continuous monitoring 
of the flux of energy and/or particles. They are rather 
based on two projective quantum measurements of the 
relevant observables (Hamiltonian, number of particles), 
performed at to and tf, where the exchanged energy and 
number of particles are obtained as differences of the out- 
comes of these two measurements @, El-El- Evidently, 
this scheme presents practical difficulties, as it is not 
clear, for example, how one could measure the macro- 
scopically large number of electrons in a lead with the 
required single electron resolution. Therefore, as other 
authors also have pointed out (see the conclusions in Ref. 
[lij]). it is necessary to bridge the gap between theory and 
experiment by extending the theory in such a way as to 
account for the possibility of continuously monitoring a 
specific quantum observable of the system. 

In this Letter we develop a multiple measurement ap- 
proach to quantum fluctuation theorems that extends the 
previous two-measurement fluctuation theorems [H, Ull - 
!13| allowing for the possibility to perform measurements 
of any observable within the interval (to,tj). Most im- 




FIG. 1: (Color online) Schematic representation of a bidirec- 
tional counting statistics experiment. Two leads (large semi- 
circles) with different electronic chemical potentials (/ii 7^2) 
and same temperature (/3i = P2) are connected through a 
double quantum dot (small circles), whose quantum state is 
continuously monitored. The state (1,0), i.e., "one electron 
in the left dot, no electrons in the right dot" is depicted. 
The transition from this state to the state (0,1) signals the 
exchange of one electron from subsystem 1 to subsystem 2. 
Hi t 2 and A/1,2 denote the Hamiltonian and number of elec- 
trons operators of each subsystem, respectively. 



portantly, we demonstrate how to use this approach to 
overcome the major problem of measuring total energies 
and particle numbers of large reservoirs in case of a trans- 
port problem as the one described in Ref. The 
salient point here is to consider the fluxes through an 
interface rather than absolute numbers of particles and 
amounts of energy. This provides a coherent and general 
theoretical framework within which the recently reported 
electron transport experiment of Ref. (lo| can be prop- 
erly analyzed and understood, and new experiments can 
possibly be devised. 

Multiple measurements scheme.— We consider a quan- 
tum system composed of subsystems, whose mutual in- 
teraction is turned on only within the finite time interval 
(to,tf). They are initially disconnected and separately in 
equilibrium at different temperatures and chemical po- 
tentials. On coupling them through a time-dependent 
interaction leads to energy and particle transfers. While 
this exchange of energy and particles occurs a quantum 
observable is monitored, as in the experiment reported 
in 0, see Fig. 1. 

We begin by considering only one intermediate mea- 
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surcment occurring at t\ E (to,t/)- The total Hamil- 
tonian H(t) is: H(t) = ^ i H i for t ^ (t ,tf), and 
H(t) = £\ H t + V(t) for t E (to,*/), where V(t) is a time 
dependent interaction term that couples the otherwise 
uncoupled subparts of the system within the time inter- 
val (to,tf). Due to the time dependence of the Hamilto- 
nian, the microreversibility principle does no longer apply 
in the way it does for autonomous systems. Under the 
assumption, however, that at each instant of time t the 
total Hamiltonian commutes with the antiunitary time- 
reversal operator 0, [H (t), 0] = 0, a generalized form of 
microreversibility is obtained which jointly uses the op- 
erator and the temporal inversion of the protocol H(t) 
[HI, EH (see the description of the reversed protocol 
below). 

For t < to each subsystem is assumed to be at equi- 
librium with an inverse temperature and chemical po- 
tential jUj of a particle species that may be exchanged 
between the subparts in presence of interaction. That is, 
at t < to the density matrix is given by the direct product 
of two grand canonical states: p = Yl { e-ftl-^*-/ 1 ^-**^ 
with <I>i the grand canonical free energy of subsystem i, 
and Mi the number operators of those particles that can 
be exchanged between the subsystems once the interac- 
tion is turned on. In absence of interaction the particle 
numbers in each subsystem are assumed to be conserved, 
i.e., [Hi, Mi] =0. In presence of interaction the total par- 
ticle number M = ^ Mi is still assumed to be conserved, 
i.e., [H(t),M] = 0. 

At time to a joint quantum measurement of all Hi's and 
Mi's is performed. As a consequence the total system's 
wave function collapses onto a common eigenstate \ipn) 
of all these commuting observables, whose eigenvalues 
are given by: Hi\ip n ) = E^\tp n ) and Mi\i/} n ) = N^ n ). 
The system then evolves according to the unitary time 
evolution Ut 1 .t governed by H(t) until time t\ when an 
observable A is measured. We assume for simplicity that 
the eigenvalues a r of A are nondegenerate and that A 
commutes with the time reversal operator 0. As a con- 
sequence of this measurement, the system's wave func- 
tion instantaneously collapses onto an eigenstate \a r ) of 
A. The evolution generated by H(t) then continues until 
the time tf, when a measurement of all Hi's and Mi's 
is again performed, collapsing the wave function to a 
joint eigenstate \ip m ) of Hi and Mi with eigenvalues E l m 
and N^ a , respectively. Accordingly, the changes of energy 
AE, = E l rn -E l n and of particle numbers AN t =N l m -N l n 
are recorded. During each realization of this protocol a 
three-state transition of the type \tpn) — > \a r ) — > |VVra) is 
recorded. This three-state transition occurs with proba- 
bility: 

P[n,r,m] = p(m, tf\r, h)p(r, ti \n, t )p nn (1) 

where 

Pun = (iMplVn) = \{e~^-^-^ (2) 



is the probability to find the system in the state \ip n ) 
at to, and p(r, t\\n,to) (p(m,tf\r,ti)) is the conditional 
probability to find the system in \a r ) at ti (|V>ro) at tf) 
provided that it was in \ip n ) at to {\a r } at ti), i.e.: 

p{m,t } \r,t 1 ) = \(^ m \U tf , tl \ar)\ 2 (3) 
p(r,*iMo) = Ka r |£/ tl , t0 |?M| 2 (4) 

Accordingly, the joint forward (F) pdf of energy and par- 
ticle number exchanges with one interruption at t\ be- 
comes: 

Pj? ({AEA, {AN,}) = E II S ^ E * E rn + K) 

m,n i 

xSiAN, - N l m + N l n )p^{m,t f \n,to)p nn (5) 

where, 

p'^m^/ln,^) = y^p(m, tf\r, ti)p(r, *i|n, t ) (6) 

r 

denotes the conditional probability of finding the state 
\ipm) at tf provided that the state \ip n ) was found at to, 
and the observable A was measured at t\. The sum in 
([6|) runs over all eigenstates of A, and the superscript t\ 
in Eqs. (|5|6p indicates the measurement occurring at t\. 

We now consider the reversed protocol specified by 
preparing the system in the state p introduced above, 
changing the Hamiltonian in time according to the re- 
versed schedule H(t) = H(tf +to — i), and performing 
measurements of a) all H^s and Mi's, at time to, b) the 
observable A at ti = to + tf — ti, c) all Hi's, Mi's, at 
time t / . During each realization of this backward proto- 
col a three-state transition of the type Q\ip m ) —> 0|a r ) — > 
Q\ip n ) is recorded with probability: 

P[m,r,n] = p(n,t f \r,t 1 )p(r,t 1 \m,to)p mm (7) 

where p(r,ti\m,t ) = \{a r \Q^U Ti to <d\ip m )\ 2 , and 

p(n,t f \r,ti) = |('0ri|6 t C/ t/ ^ 1 6|ar)| 2 , with U t > , t the time 

evolution operator governed by H(t). Thus the backward 
(B) pdf of energy and particle number exchanges, with 
an interruption at t\ is: 

Pf ({AE,}, {ANi}) = E II S ( AE * E n + E l) 

n,m i 

xS{AN t -N l n + N'jf 1 (n,t f \m, t )p mm (8) 

where 

V 

is the conditional probability to find the system in the 
state \ip n ) at time tf provided that it was in \ip m ) at to 
and the observable A was measured at t\. 

Using [H (t), 0] = 0, and expressing the time evolution 
operator as a time ordered product one finds, in a similar 
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way as in that 0tZ7 F = U tlttj , and 0t[/ t ~ 

UtoM- Thus p(n,tf\r,t!) = p(r,ti\n,t ), p(r,ti\m,t ) 
p(m,tf\r,ti), and consequently 



p 1 (m,tf\n,t ) = p 1 (n,tf\m,t ) 
follows. Combining Eqs. (|2I5I8I10[) we obtain: 



if ({-ASi},{-AJVi}) 



(10) 



(11) 



which reads exactly as the two-measurements fluctua- 
tion theorem [12], with the major difference that now, 
a third projective measurement is performed at t\. Here 
the free energy difference does not appear in Eq. (fTTj) 
because we assumed if (to) = H(tf). We stress that 
the backward and forward probabilities with intermedi- 
ate measurement Pp , Pp , in general differ from the cor- 
responding probabilities without intermediate measure- 
ment Pp, Pb of Ref. 12J. Remarkably, however, their 

ratio remains unaltered Pp 1 /Pp = Pp / Pb- In both cases 
though, particle numbers and energy measurements are 
assumed to be performed at to and t/. As discussed in 
the introduction, these measurements are practically im- 
possible when the subsystems are macroscopic reservoirs, 
as is the case with quantum transport problems. In the 
following we will show a way to get around this problem 
by taking advantage of the possibility opened by Eq. (fTTj) 
of performing intermediate measurements. 

The salient point in the derivation of Eq. (fTTj) is 
that the generalized time reversal invariancc property 
of the transition probability without interruptions, i.e., 
p(m, tf\n, to) = p(n,tf\m,to) [l2j], continues to hold, Eq. 
(|10[). even if a projective measurement is performed at t%. 
Evidently the result (fTTj) can be extended to the case of 
multiple intermediate quantum measurements of possibly 
distinct observables A k , k = 1 . . . K occurring at times t k 
during the forward protocol, and at times tk = to +t/ — tfc 
during the backward protocol. Thus we conclude that the 
fluctuation theorem is unaffected by the action of inter- 
mediate projective quantum measurements. 

Indeed a more detailed result holds. Consider the 
joint probability for a sequence of outcomes no,rTk,nf = 
no, Tii ■ ■ -Tin,Tif, and the joint probability of the reversed 
sequence n/, n^, no = nf,n n ...ni,no for the backward 
protocol: 

K 



P[n ,n k ,n f ] = p(n k +i, t k +i\rik, tk)pn n (12) 

fe=0 
K 

P[nf,rTk\,n Q } = JJ p(n k , tk\n k +i, tk+i)p nf n f (13) 



fc=0 



where tf t o — to./, and no,n/ label common eigenstatcs 
of {Hi, H 2 ,Afi,Af2}, and n k labels eigenstates of Ak. 
Using [H(t),Q] = one finds QUrr fc 0t = U tk+utk , 



which leads to p{n k , t k \n k+ i,t k+ i) = p(n k+ i,t k+ i\n k ,t k ). 
Then, using Eqs. (|2I12I13|) . we readily find the following 
detailed fluctuation theorem: 

P[n,n k ,m] 



P[m,nfc, 



n 



A(- E m-- E »-« JV m+/- l i iV ») 



(14) 



Notably, the right hand side depends only on initial and 
final eigenvalues of Hi, Mi, but not on the eigenvalues of 
Ak- Eq. (|14[) can be regarded as the quantum mechan- 
ical version of a formula that was first put forward by 
Bochkov and Kuzovlev for classical trajectories (see Eq. 
(7) of Ref. 0|; a related result was derived recently for a 
classical master equation in [Hj]). 

Application to counting statistics experiments.— In the 
counting statistics experiment of Ref. [TcJ two leads hav- 
ing the same temperature j3 are connected via a double 
quantum dot. We identify the left reservoir plus left dot 
as subsystem 1, and the right reservoir plus right dot as 
subsystem 2, see Fig 1. The two subsystems are coupled 
through a small (but important) timc-indcpcndcnt inter- 
action term V. At time to an electric potential difference 
Acp is applied across the two subsystems. We model this 
by saying that immediately after the time to the density 
matrix is the factorized grand canonical state p speci- 
fied above with eAip = fj,i — /i2- Afterwards the system 
follows the evolution according to H = Hi + H2 + V. 
This is equivalent to assuming that V is switched on in- 
stantaneously at to. Since V is small the energy of the 
total system made of the leads and the double quantum 
dot hardly changes, hence, AE\ ~ — AE%. Due to the 
conservation of the total number of electrons, we have 
ANi = —AN 2 = q. A time-reversal invariant observable 
A, detecting whether there is an excess electron in either 
of the two dots is continuously monitored, resulting in 
a signal nt. The total number of exchanged particles is 
calculated from vX, by means of a function q[nk] that de- 
termines the number of transitions from the state (1,0) 
(i.e., one electron in the left dot and no electrons in the 
right dot) to the state (0, 1) (i.e., no electrons in the left 
dot and one electron in the right dot), minus the num- 
ber of opposite transitions. Evidently q[nt] = — q[nk] 
holds. The only difference with the multiple measure- 
ment scheme presented above is that in the experiment 
no measurements of total particle numbers in the two 
subsystems are performed at to and tf. However since 
the initial density matrix is a mixture of states with defi- 
nite particle numbers in each subsystems, a measurement 
at to can be abandoned if one is only interested in rel- 
ative changes of particle numbers. By means of regis- 
tration of each transition between the two subsystems, 
not only the flux of electrons is monitored but it is also 
guaranteed that during the whole process the particle 
numbers of each subsystem have definite values. There- 
fore also at tf the state of the total system is a mixture 
of states with definite particle numbers, hence a mea- 
surement of particle numbers need not be performed at 
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tf. Thus g[n£] coincides with the value that would have 
been obtained from electron numbers measurements, i.e., 
= N} n -N^ = -N^ + N%. Therefore the pdf of 
particle exchange is: 



P 



(?) 



6(q-q[nt\)P[ni 



(15) 



where we introduced the marginal probability P[nt} = 
^ m P[n,n||m]. Since the backward and forward pro- 
tocols coincide in this case, we havePj tfe} = P|* fe} . Using 
Eqs. (|14I15|) with q[nt] = —q\nk], wc conclude that the 
pdf of electron exchange indeed satisfies: 



,{**} 



(<?) 



(16) 



in good qualitative agreement for small small r, and ex 
cellent quantitative agreement for large r, with experi 



mental findings of Rcf. [10(, see Fig. 2a therein [17|. 
The same experimental setup could be used to confirm 
the validity of the detailed fluctuation relation 



PH] = P[fck]e 



P(fn-fi2)q[n k ] 



(17) 



that follows from (|T4j) . 

Remarks.— One assumption of our theory is that the 
initial state of the system is the direct product of two 
grand-canonical states. For this particular choice of the 
initial state the fluctuation theorems (|16I1T[) do hold in- 
dependently of the duration r of the protocol. In exper- 
iments, however, unavoidable interactions may lead to 
correlations between the sub-parts of the system already 
before the protocol has started. These initial correla- 
tions in general will lead to transient deviations from the 
fluctuation theorems which then arc expected to hold as 
steady state fluctuation theorems in the limit of large r, 
as is the case with experiment fioj . 

The assumptions that the measured observables 
have non-degenerate eigenvalues can be abandoned with- 
out any change of our central result formulated in Eqs. 
(fTTI H"4")) . We further found that our main result remains 
true if the quantum measurements are performed with re- 
spect to positive operator valued measures (POVM's) in- 
stead of von Neumann projection valued measures. This 
is not so surprising after all, being the former less invasive 
than the latter. The main result holds for any number 
of subsystems. In particular for a single closed system it 
gen eralizes the Tasaki-Crooks work fluctuation theorem 



20l | . In case of two susbsystems with equal temper- 
ature and no matter exchange it generalizes the Tasaki- 
Crooks fluctuation theorem for open systems [2l|, l22j ]. 
Multiple measurements of reservoir energies (not fluxes) 
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were considered for this case in Rcf. 

Far from being a problem of purely academic interest, 
the multiple measurement scheme pursued here, helps 
bridging the existing gap between theory and experi- 
ments concerning quantum fluctuation theorems. The 



measurement of energy and particle content of the sub- 
parts of a system may be practically impossible, espe- 
cially when the subparts are macroscopic objects like 
the leads in an electron counting statistics experiment. 
On the other hand energy and matter exchanges may be 
measured by monitoring the respective fluxes through in- 
terfaces, making the cumbersome measurements of total 
energies and particle numbers obsolete. Here we have 
shown that the fluctuation theorem continues to hold in 
this case of continuous monitoring, Eqs. (|11I14|) . pro- 
vided a theoretical explanation for an experimentally ob- 
served relation, Eq. (jT6|) , and predicted that a more de- 
tailed relation, Eq. (fl7j) should hold in the same set-up. 
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